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Abstract 

The family of complex PT-symmetric sextic potentials is studied to show that for 
various cases the system is essentially quasi-solvable and possesses real, discrete energy 
eigenvalues. For a particular choice of parameters, we find that under super symmetric 
transformations the underlying potential picks up a reflectionless part. 
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1 Introduction 

Searching for non-Hermitian PT-symmetric Hamiltonians has acquired much interest in 
recent times (see e.g. fl], |2|, |3], |J and references quoted therein). For one thing, a rather 
large subclass of such Hamiltonians has been found to possess real energy eigenvalues. For 
another, in at least some cases, it is seen that a complex shift of the coordinate x G (— oo, oo) 
does not affect the overall normalizability of the wave functions, while at the same time 
retaining the real character of the energy spectrum. 
The purpose of this letter is twofold: 

(i) We examine the general problem of a complex sextic potential from the point of view of 
determining exactly a finite number of eigenvalues and eigenfunctions. A suitable ansatz 
scheme leads us to find discrete real energy levels under quite general conditions. 

(ii) We point out that some of our results can also be arrived at by performing a complex 
shift of the coordinate on the reduced sextic potential consisting of even-power terms only. 
However, our results cover a much greater ground. In particular, we find it possible to 
generate an additional complex reflectionless term in the sextic potential by employing 
supersymmetric transformations. 

2 Complex sextic potentials and their solutions 

To get started, let us consider the following general representation of a sixth-degree potential 

V{x)=j2c lX \ (1) 
1=1 

satisfying the Schrodinger equation (in units h = m = 1) 




(2) 



where for V(x) to be PT symmetric, Ci, C3, C5 G iK, but C2, C4, Cq G E. 



We make the ansatz that the wave function is of the form 




(3) 



2 



where f(x) is some polynomial function of x, which, for complex potentials, is typically of 
the type J2m=o a m{^) m . For the real analogue of ([I]) consisting of even power terms only, 
f(x) is known to have a given parity. 
We focus on the following choices of /: 

(a) f(x) = l, 

(b) f(x)=x + a , 

(c) f(x) = x 2 + a\X + a , 

but can generalize to higher degrees as well. For complex potentials, clq is imaginary in (b), 
whereas a± is imaginary, but a is real in (c). 

Without going into the details of calculations, which are quite straightforward, let us 
summarize our results. 

2.1 The f{x) = 1 case 

Here the potential parameters are found to be related to the fe's as 

ci = -36 3 + 26i6 2 , c 2 = -664 + 36i6 3 + 26;;, c 3 = 46i6 4 + 66263, 

c 4 = 86264 + pl 05 = 126364, c 6 = 8bl (4) 

Without loss of generality, we can choose c 6 = | to fix the leading coefficient of V(x). 
It gives 64 = ±|. We take the positive sign to ensure normalizibility of the wave function, 
which reads 

ip(x) = exp (— bix — b 2 x 2 — 6 3 x 3 — \x^ . (5) 
The associated energy level is given by 

E = 62 - \b\. (6) 

Now 61 and 63 imaginary make c±, C3, C5 imaginary too, so we have a complex PT- 
symmetric potential. Note that the energy eigenvalue in such a case is real, and the corre- 
sponding wave function is PT-symmetric. 
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2.2 The f(x) = x + a case 

The wave function is of the form 

ip{x) = (x + ao) exp (— b\x — b2X 2 — b 3 x 3 — |x 4 ) (7) 

for 

ci = -663 + 2&!& 2 + a , c 2 = -f + 36i6 3 + 2&2, c 3 = 61 + 66 2 &3, 

C4 = 26 2 + |62, c 5 = 36 3 , ce = i (8) 

There is also a condition on ao, 

a 3 - 3b 3 a 2 + 2b 2 a - h = 0. (9) 

The energy is given by 

E = -\b\ + 36 2 - 3a 6 3 + a\. (10) 
Let us discuss some important special cases of this scheme. 

2.2.1 6x = 63 = 

The condition (g) reduces to 

a (a 2 + 26 2 ) = 0. (11) 

(i) If ao = 0, then there is no imaginary term in the potential, and this corresponds to 
the n = 0, negative-parity result of ref. |J. The energy eigenvalue is given by E = 3b 2 and 
shows a single level. 

(ii) The other solution of flTT|), namely = — 26 2 can be studied according to whether 
al > or ag < 0. 



If ag > 0, then a linear term is present in V(x) with c\ = ±y2|6 2 |, c 2 = 2b 2 . — ~, and 
C4 = 26 2 - Of course C3 = C5 = 0. The energy is E = b 2 < 0. Thus we have two different 
real potentials with the same energy eigenvalue. The linear term breaks PT invariance of 
the potential and the wave function as well. So, in this respect, ao real can be viewed as 
an explicit symmetry breaking parameter. 
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On the other hand, if < 0, we get two different complex potentials, corresponding to 



c\ = ±iy/2b 2 , with the same real energy eigenvalue: 



V(x) = l x s + 2b 2 x 4 + (2b 2 2 -l) x 2 ±i x /2b 2 x, 



E 

ip(x) 



b 2 >0, 



x ± iJ2b 2 ) exp (— b 2 x l — \x A 



(12) 
(13) 
(14) 



The potential ( |T2"D is PT symmetric, while the wave function (|T^) is odd under PT symme- 
try. 

2.2.2 b x = 0, b 3 ^ 

The solution for ao = turns out to give the same conclusions as previously obtained. The 
second solution 

a = \ (363 ± ^961 - 86 2 ) (15) 

yields two possibilities according as 63 G K or 63 G iR. If 63 G R, we have b\ > ^b 2 , implying 
two different real potentials with the same energy eigenvalue, except for the equality sign 
in (HH). 

If however 63 G iR, then ao must be imaginary with b\ = — \b 3 \ 2 < |6 2 . Here too we have 
two possibilities of obtaining two different complex potentials with the same real energy 
eigenvalue, except for the equality sign in (|TJ). 



2.3 The f(x) = x 2 + cl\X + a case 

The complete set of solutions leading to more than one energy level corresponds to 

1 



and gives 
V(x) = 

E± = 

v>± = 



ai = 26 3 



a 



2b 2 -b'i± V(26 2 - 3&§) 2 + 2 



|x 6 + 36 3 x 5 + (2b 2 + pi) x A + 26 3 (4& 2 - bj) x 3 + [2 (b 2 2 + 3b 2 bj - 3bj) - 



+ b 3 (4b 2 2 -4b 2 b 2 3 -7)x 



2bUb 2 -bl) +36 2 -&3± V(2&2-36 



x 2 + 2b 3 x + \ [2b 2 - b 2 3 T V ( 2 &2 - 3&§) 2 + 2 
x exp [-26 3 (b 2 -bf)x - b 2 x 2 - b 3 x 3 - \x A 



(16) 

x 2 
(17) 
(18) 

(19) 



The results ([T6|) — (|Tg|) are valid both for real and imaginary fe's. In the latter case, PT 



symmetry is good for the potential and the wave function, whereas in the former one has 
a symmetry breaking. Note that in both cases E + > E_. 

Concerning the case where 63 is imaginary, we have a complex PT-symmetric two- 
parameter family of potentials with two distinct real energy levels. This is truly a non- 
trivial result and puts the spirit of quasi-solvability in the complex domain. Indeed, for 
the particular case of 63 = and 62 = 7/2, we recover the n — 1, positive-parity results 
of the one-dimensional even-power potential V(x) = \x % + 7X 4 + | (7 2 + /i) x 2 , where /i = 
—3 — An — 2r and r is associated with the (— l) r parity of n + 1 levels. 

The converse also works. The results (|16D-(|19|) can be derived from the even-power 
sextic potential of ref. by a translation x — > x + 63. If 63 is imaginary, then this 
translation amounts to a complex shift, whose viability has already been pointed out in 



In order to enlarge the class of PT-invariant potentials, SUSY methods have been used 
thoroughly [Q, |J]. Here, we also outline the procedure to generate superpartners of the po- 
tentials considered in the previous section, which share their PT properties. The procedure 
is based on the construction of a superpotential W(x) = —ip'(x)/ip(x), and therefore is 
sound in so far as the logarithmic derivative of ip(x) is well behaved on the real axis. 

One can start from the wave function ([5]), the parameters bi and 63 being assumed 
imaginary. Then the superpartner 



refs. 



3 A supersymmetric viewpoint 




(20) 



becomes again a sextic potential, and there is no real enlargement. 

More interesting is the case where one starts from the wave function (fb3| ) and constructs 
the partner of (0). By substituting ([TJ]) into (PUD, we find V to be 




(21) 
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where the piece in parentheses is clearly reflectionless. The latter is indeed reminiscent of 
the "transparent" complex potential obtained in ref. |2[], which is invariant under PT and 
gives a trivial S-matrix. Moreover this piece has an associated zero-energy bound state, 
given by ^o(x) = C (x ± iy/Tj^j , where C is a constant. [] 

The form ([H]) also generalizes the result obtained for the harmonic oscillator in ref. |3| 
to the sextic case. The fact that 62 7^ makes the potential rather appealing in that if 6 2 
were vanishing, there would be a singularity on the real axis and the potential would be 
rendered ill defined. Here, because of a shift of the singularity, Q2T| ) remains well posed in 
the complex plane cut from x = ^\\f%~2 to x — =Fioo, respectively. 

Finally, one can consider the wave function ([19]) and construct the partner of (Jl7p . 



Again one can see that V(x) is not trivial: for 63 imaginary, ip'(x)/i/j(x) is well behaved on 
the real axis. 

All results we have obtained can be generalized to the case where one shifts the variable 
x by a translation b, with b real. In such a case, the potentials and wave functions are 
reparametrized correspondingly. One may worry however about PT properties since one 
can generate subleading powers in x from a given power. However, one should realize that 
now the parity operation can be defined with respect to a mirror placed at x — —b, so that 
x + b = x — (—6) goes to — x — 2b— (—6) = —(x + b). When this translation b is performed, 
the reflectionless potential contained in fl2ll) becomes precisely that considered in ref. ||. 

4 Conclusion 

To conclude, we have solved a complex PT-symmetric sextic potential in its most general 
form within a suitable ansatz scheme for the wave functions and shown how the associated 
energy levels turn out to be real. We have also demonstrated using SUSY the possibility of 
generating a complex reflectionless part in the potential. Although we have restricted our 
discussion up to the quadratic order in the coefficient of the exponential representing the 
wave function, it is obvious that we can build up, in an identical way, higher-order states. 



lr This transparent potential may also be viewed as the particular case I — 1 of the generic potential 
7^1(1 + 1) (x ± i\/2&2) that is typical of a centrifugal barrier in a radial context (but with no singularity) 
and which has a zero-energy bound state ^o(x) = C (x± 1^/282) ■ 
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